Configurational Temperature of Charge-Stabilized Colloidal Monolayers 
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Recent theoretical advances show that the temperature of a system in equilibrium can be mea- 
sured from static snapshots of its constituents' instantaneous configurations, without regard to their 
dynamics. We report the first measurements of the configurational temperature in an experimental 
system. In particular, we introduce a hierarchy of hyperconfigurational temperature definitions, 
which we use to analyze monolayers of charge-stabilized colloidal spheres. Equality of the hyper- 
configurational and bulk thermodynamic temperatures provides previously lacking thermodynamic 
self-consistency checks for the measured colloidal pair potentials, and thereby casts new light on 
anomalous like-charge colloidal attractions induced by geometric confinement. 



The temperature of an equilibrium ensemble of parti- 
cles is denned conventionally in terms of the particles' 
mean kinetic energy, without regard for their instanta- 
neous positions. In 1997, Rugh pointed out that the 
temperature can also emerge from other ensemble aver- 
ages over geometrical and dynamical quantities 0] . This 
notion is expressed more generally 0, as 
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where angle brackets indicate an ensemble (or time) av- 
erage, T = {qx, . . . , q 3 N,pi, ■ ■ ■ ,P3n} is the instanta- 
neous set of 3N generalized coordinates qj and their 
conjugate momenta pj for an TV-particle system, TL = 
E^=i Pj I (2m.) + V({qj}) is the Hamiltonian associated 
with the conservative JV-particle potential V({qj}), and 
B(r) is an arbitrary vector field selected so that both the 
numerator and denominator of Eq. Q are finite and the 
numerator grows more slowly than in the thermo- 
dynamic limit. Choosing B(T) = {0, . . . ,0,pi, ■ ■ -P3n} 
yields the usual kinetic definition of temperature. Choos- 
ing instead B(T) = — W({qi}) yields a formally equiva- 
lent result, 
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which depends only on the particles' instantaneous con- 
figuration, and not on their momenta. Variants of Eq. (J2J 
have been widely adopted in molecular dynamics simu- 
lations 0, Q, but have not previously been applied to 
experiments. 

In this Letter, we use the configurational temperature 
formalism to probe macroionic interactions in monolayers 
of charged colloidal spheres dispersed in water and con- 
fined between parallel glass plates. These measurements 
provide sensitive self-consistency tests for the measured 
inter-particle pair potentials, and thereby provide new in- 
sights into the long-standing conundrum of anomalous at- 
tractions in geometrically confined charge-stabilized dis- 
persions. 



Directly applying Eq. J5J requires the full iV-particle 
free energy, which is rarely available. Simplified forms 
emerge for systems satisfying certain conditions. For ex- 
ample, if V({qi}) is the linear superposition of pair po- 
tentials, it(r), then Eq. J2J reduces to |4(, 
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where Fj = — Ej^i ^i u ( r ij) is the total force on particle 
i due to its interactions with other particles, is the 
gradient with respect to the i-th particle's position, r^, 
and Tij — |r,; — | is the center-to-center separation be- 
tween particles i and j. The temperature is reflected in 
the instantaneous distribution of forces because objects 
explore more of their potential energy landscape as the 
temperature increases. 

Equation J2J) may be generalized into a hierarchy of 
hyperconfigurational temperatures by choosing B(r) = 
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for s > 0. These higher moments are more sensitive 
to the input potential's detailed structure than T con F = 
(2) . They also can be applied to three-dimensional sys- 

iK is 



terns with long-ranged 1/r potentials, for which T c , 
ill-defined. 

The same hierarchy also may be derived from the clas- 
sical hypervirial theorem ({/, Ti}) = 0, where {• ■ ■ } 
is the Poisson bracket, by selecting f(T) = Ej=iPi Fj- 
More generally, f(T) can be any finite- valued function 
that does not explicitly depend on time. 

For systems with short-ranged potentials, dropping ad- 
ditional terms of 0(1 /TV) from Eq. {IJ yields @]: 
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the second of which is proposed here for the first time. 
These definitions' different dependences on sample size 
N are useful for comparison with TJ" . 

Temperature definitions based on configurational in- 
formation are ideal for studying colloidal spheres dis- 
persed in viscous fluids such as water. The fluid acts 
as a heat bath at temperature T. It also randomizes 
the particles' motions over intervals longer than the vis- 
cous relaxation time, typically measured in microseconds. 
As a result, the colloids' instantaneous momenta are not 
easily accessible. Their positions, however, are readily 
measured using standard methods of digital video mi- 
croscopy 0. Calculating the temperature from the re- 
sulting positional data then requires accurate knowledge 
of the colloids' interactions. 

The mean-field theory for macroionic interactions Q 
predicts that charge-stabilized colloidal spheres should 
repel each other through a screened-Coulomb potential: 
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where /9 _1 = ksT is the heat bath's thermal energy 
scale, a is the spheres' diameter, r is their center-to- 
center separation, Z is the effective charge number on 
each sphere and Xb — e 2 / (A-KeksT) is the Bjerrum length 
in a medium of dielectric constant e. The Debye-Hiickel 
screening length k _1 sets the range of the effective elec- 
trostatic interaction and depends on the concentration c 
of (monovalent) ions through k 2 = 4itcXb- 

Despite its success at explaining bulk colloidalphe- 
nomena mean-field theory qualitatively fails to 
explain the strong and long-ranged attractions observed 
when charged spheres are confined between parallel glass 
walls [a, llCl 111! |l2( . The crossover from monotonic re- 
pulsion to long-ranged attraction with increasing confine- 
ment is demonstrated in Fig. Efforts to explain this 
anomaly through non-mean-field mechanisms so far have 
not yielded the experimentally observed effect [13j . Given 
the apparent difficulty of explaining anomalous attrac- 
tions on the basis of colloidal electrostatics or electrody- 
namics, various other explanations have been proposed, 
most of which focus on possible experimental artifacts. 
For example, nonequilibrium hydrodynamic coupling has 
been shown |14j to explain one measurement based on op- 
tical tweezer manipulation [l5j . but cannot be relevant 
for measurements based on video microscopy of colloid 
in equilibrium 0,0, Concern also has been 

raised that such imaging measurements can fall victim 
to correlated artifacts due to statistical fluctuations 0] • 
We have shown, however, that anomalous attractions still 
are clearly resolved given adequate statistics . A still 
greater concern is that uncorrected many-body artifacts 
could mimic attractions in purely repulsive dispersions. 
Such artifacts certainly can arise at hi gh d ensities, as 
has been demonstrated experimentally |l8| . Evidence 
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FIG. 1: Pair potentials u(r) for silica colloid a = 1.58 /im 
in diameter extracted from the measured radial distribution 
functions g(r) using the HNC approximation for suspensions 
at H = 9 fim, na 2 — 0.0654 (squares) and H = 195 /im, 
na 2 = 0.0797 (circles). Solid curves are, respectively, fits to 
a 5-th order polynomial and to Eq. Q (k _1 = 180 ± 10 nm, 
Z = 6500 ± 1000). 



for pairwise additivity at lower densities has relied prin- 
cipally on comparisons over a range of concentrations, 
for which variations in chemical environment could mask 
other effects. 

We resolve all such ambiguities by exploiting the con- 
figurational temperatures' sensitivity to inaccuracies in 
the input potential as self-consistency tests for colloidal 
interaction measurements. In particular, we use pair po- 
tentials measured according to Refs. [g, ^3 an d 01 
to compute the configurational temperature of geometri- 
cally confined colloidal monolayers and compare the re- 
sults with the bulk thermodynamic temperature, with 

( s) 

' = T signaling thermodynamic self-consistency. De- 
viations could result from a breakdown of pairwise ad- 
ditivity in a system with non-trivial many-body interac- 
tions, a departure from equilibrium in a system subjected 
to hydrodynamic forces, or simply inaccuracy in u(r). 

Our samples consist of uniform silica spheres a = 
1.58 ± 0.03 fim in diameter (Duke Scientific Lot 24169) 
dispersed in deionized water and loaded into hermetically 
sealed sample volumes formed by bonding the edges of 
glass #1.5 coverslips to the surfaces of glass microscope 
slides. The separation H between these surfaces is set 
during construction and establishes the degree of confine- 
ment. Reservoirs of mixed bed ion exchange resin main- 
tain the ionic strength below c « 10 /LtM, correspond- 
ing to a Debye-Hiickel screening length k- 1 > 150 nm. 
Under these conditions, both the silica spheres and the 
glass walls develop negative surface charge densities of 
roughly 10 3 eo /im" 2 , where eo is the electron charge 
[l9) . Samples are mounted on the stage of a Zeiss Ax- 



3 



iovert 100 STV microscope, after which the dense sil- 
ica spheres rapidly sediment into a monolayer roughly 
h = 900 nm above the coverslip with out-of-plane fluc- 
tuations smaller than 300 nm. The bright-field imaging 
system provides a magnification of 212 nm/pixel on a Hi- 
tachi TI-11A monochrome CCD camera. The resulting 
video stream is recorded and digitized into 60 deinter- 
laced video fields per second, each of which is analyzed 
to yield the instantaneous distribution 
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of N(i) spheres at locations rj(f) at time t. Measure- 
ments were performed at room temperature over peri- 
ods of roughly one hour, with temperature fluctuations 
smaller than ±1 K over the course of any measurement. 

From p(r,t), we calculate the radial distribution func- 
tion, 
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where n = N/A is the areal density of N = (N(t)) par- 
ticles in the field of view of area A, and the angle brack- 
ets indicate averages over angles and time. Assuming 
that the system is in equilibrium and that its interactions 
are isotropic and pairwise additive, g(r) may be inverted 
to obtain an estimate for u(r) [lfllll|]. At low enough 
densities, the Boltzmann distribution provides the nec- 
essary relationship: (3u{r) — — lim rw o lng(r). Unfor- 
tunately, the need to sample g(r) adequately in a lim- 
ited field of view requires higher concentrations 0, ^3 > 
with 0.05 < na 2 < 0.1 being typical for our experi- 
ments. We extract candidate pair potentials u(r) such 
as the examples in Fig. ^ using the liquid-structure in- 
version method based on the Ornstein-Zernike equation 
with hypernetted-chain (HNC) and Percus-Yevick (PY) 
closures 0, 0, . As we have reported previously, the 
scdimcntcd silica spheres repel each other as predicted by 
Eq. J7J) at the larg est inter-plated separations considered, 
H = 200 /jm [13, Uj| ■ Reducing H does not perceptibly 
change the spheres' equilibrium height h above the lower 
glass wall, yet nonetheless introduces a minimum into 
u(r) consistent with a long-ranged attraction 0. The 
example in Fig. 2] at H = 9 /im has a minimum roughly 
0.3 ksT deep at r = 1.5 a. Well-resolved minima are 
evident for spacings as large as H < 30 /im. 

Even the small amount of scatter in the measured pair 
potentials yields unacceptably large fluctuations in the 
derivatives used to calculate the configurational temper- 
atures. We avoid this by fitting the experimental data 
to fifth-order polynomials, as shown in Fig. ^ and using 
the fits as inputs to Eqs. |@J, © and 

In calculating the total force Pi on the i-th. particle, 
we must consider all relevant pair interactions. Particles 
close to the edge of the field of view may have strongly 
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FIG. 2: (a) Estimating the interaction range for the H = 9 fim 
data set using T(r). Configurational temperatures were cal- 
culated over 0.93 < r/a < 2. (b) Configurational and hy- 
perconfigurational temperatures as a function of sample size, 
showing finite-size scaling through fits to second-order poly- 
nomials in 1/N (third-order for T con 2). 



interacting neighbors just out of view. To avoid errors 
due to the resulting spuriously unbalanced forces, we 
restrict ensemble averages to particles farther from the 
edges than the range R of the interaction. We estimate 
R by plotting T(r)/T = 2n(r / a) g(r)\Vu{r)\ 2 /V 2 u{r), 
as shown in Fig.|3a), which qualitatively gauges contri- 
butions to the configurational temperature due to parti- 
cles at separation r. This is most useful for systems with 
weak three-body correlations. Typically, R < 2.0 a for 
our samples. 

Sample imperfections such as a small number of dimers 
also distort the apparent force distribution We omit such 
particles from the ensemble averages. By contrast, we 
include the small number of pair interactions with r/a < 
1, which arise from the samples' polydispersity and also 
because of projection errors due to out-of-plane motions 
for particles near contact. Omitting these has little effect 
on T CO nF, but leads to large systematic errors in higher- 
order T^ s) . 

The restricted sample includes only N — 0(100) par- 
ticles, hardly the thermodynamic limit. To correct for 
finite-size effects, we deliberately sub-sample our data, 
and plot the configurational temperature as a function of 
1/N, as shown in Fig.[2Ib). Polynomial fits account for 
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FIG. 3: Configurational temperatures at different plate sepa- 
rations. The uncertainty of each extrapolated temperature is 
less than 0.02. 



the finite-size scaling of Eq. and permit extrapola- 
tions to large N. As expected, T conl and T con2 are more 

(s) 

sensitive to sample size than because their deriva- 
tions involve ignoring more terms of 0(1/ N). Despite 

( s) 

these differences, all orders of , as well as T conl and 
7con2i extrapolate to the thermodynamic temperature in 
the large- N limit. This successful outcome strongly sug- 
gests that the experimentally determined potential, in- 
cluding its long-ranged attraction, accurately describes 
the dispersion's equilibrium pair interactions. 

The configurational temperatures are exceedingly sen- 
sitive to small variations in /3it(r), with systematic ad- 
justments as small as 0.01 in the short-range repulsive 
core leading to variations in the configurational temper- 
ature as large as ten percent. Truncating the attractive 
part of u(r) increases the configurational temperatures 
by fifty percent. Similarly, failure to precisely correct 
the imaging system's aspect ratio can lead to large devi- 
ations. This potential source of error can be monitored 
by factoring T con p into components along and transverse 
to the video scan lines and comparing the results. De- 
liberately rescaling one axis by 0.01 causes the apparent 
temperatures along the two directions to differ systemat- 
ically by as much as ten percent for all samples. 

Figure demonstrates that the configurational tem- 
peratures are consistent with the thermodynamic tem- 
perature over the entire range of wall separations from 
H = 3.2 fim to H = 195 /im, despite variations in the 
form of the associated pair potential This allows 

us to draw several conclusions regarding the nature of 
confinement-mediated colloidal interactions. Primarily, 
we conclude that the measured pair potentials accurately 
and self-consistently describe the colloidal particles' in- 
teractions. To the extent that the configurational tem- 
peratures are sensitive to local departures from equilib- 
rium [13, the result T conf /T = 1.00 ± 0.04 for all of our 
samples suggests that anomalous confinement-induced 
like-charge colloidal attractions cannot be ascribed to 
nonequilibrium mechanisms such as hydrodynamic cou- 
pling due to transient flows in the sample 



In summary, we have calculated the configurational 
temperature for experimentally determined distributions 
of colloidal silica spheres using their measured pair po- 
tentials as inputs. The success of this procedure provides 
a thermodynamic self-consistency test for the measured 
potentials, and thus establishes that confinement induces 
equilibrium pair attractions between charged colloidal 
spheres. It furthermore demonstrates the configurational 
temperature to be a powerful new tool for experimental 
condensed matter physics. 

We are grateful to Owen Jepps for introducing us 
to the notion of configurational temperatures, to Sven 
Behrens for extensive discussions, and to the donors of 
the Petroleum Research Fund of the American Chemical 
Society for support. 
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